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High parton density effects with energy obey non-linear QCD evolution equations 
for which exact solutions are not known. The mathematical class to which the non- 
linear Balitsky-Kovchegov equation belongs is identified, proving the existence of 
asymptotic in energy traveling wave solutions which are "universal" i.e. indepen- 
dent of the initial conditions and of the precise form of the non-linearities. This has 
an direct impact on geometrical scaling and the diffusive transition to saturation, 
which is shown to be "normal" for constant QCD coupling and "abnormal" for 
running coupling. 

1 Introduction 

Considering the scattering of a hard projectile (e.g. a massive QCD dipole) on an 
extended target, the Balitsky Fadin Kuraev Lipatov (BFKL) £Q evolution equation 
implies a densification of gluons and sea quarks with incident energy, while they 
keep in average the same size. It is thus natural to expect a modification of the 
evolution equation towards a saturation regime. Recently, a theoretical appoach to 
saturation has been found |.'SI 1| related to non-linear evolution equations of the gluon 
density in the framework of perturbative QCD. In the transition to saturation, the 
exponential growth regime related to the BFKL kernel gets modified by non-linear 
terms, leading to the Balitsky-Kovchegov (BK) equation 0]. A more general non- 
linear functional equation is expected to take into account the multiple correlations 
and to describe the fully saturated phase [3]. The aim of our approach |5I6I7| is 
to explore the mathematical properties of the BK equation and derive its physical 
consequences for saturation in QCD. 

2 Saturation and Non-Linear Equations 

The Balitsky-Kovchegov (BK) equation 4 considers the energy evolution within 
the QCD dipole Hilbert space [8]. To be specific let us consider N(Y, xoi), the 
dipole forward scattering amplitude and define 



Jo xoi 

Within suitable approximations (large N c , summation of fan diagrams, spatial 
homogeneity), this quantity obeys (see the second reference in 0]) the nonlinear 
evolution equation 




(1) 



dyAf = a\ (-d L )J\f - aAf 2 , 



(2) 
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Figure 1. Typical traveling wave solution. The function u(t, x) is represented for three different 
times. The wave front connecting the regions u = 1 and u = travels from the left to the right 
as t increases. That illustrates how the "strong absorption" or saturated phase region invades the 
"transparency" region. 



where a — a s N c /ir, = 2-0(1) — "0(7) — "0(1 — 7) is the characteristic function 
of the BFKL kernel [Tj, L = log(fc 2 ). In a first stage [S], let us consider the kernel 
expanded up to second order around 7 = § . 

Eq. (0) boils down to a parabolic nonlinear partial derivative equation: 

dYJV=a{ X (±) + ±x"(i)(d L + i) 2 }M-aJV 2 . (3) 

The mathematical point |S] of our recent approach is to remark that the struc- 
ture of Eq.© is identical (by a suitable linear redefinition N{L,Y) — > u(x,t) and 
for fixed a) to the Fisher and Kolmogorov-Petrovsky-Piscounov (F-KPP) equation 

M 

d t u(t, x) = d 2 x u(t, x) + u(t, x)(l - u(t, x)) (4) 
which appeared in the problem of gene diffusion and annihilation (1938). 

3 Universality and Traveling Wave Solutions 

The remarkable mathematical property of the F-KPP equation is the existence of 
traveling wave solutions of the F-KPP equation '9 at large times. This means that 
there exists a function of one variable w such that 

u(t^+oo,x) ~ w{x- It- § log* + 0(1)} (5) 

uniformly in x. Such a solution is depicted on Fig.l|T|). 

This analysis can be extended [Hj to the study of the equation with the full 
kernel. Indeed, only the second-order expansion around a given critical value 7 = 
7 C = .6275... is relevant. Let us describe its general consequences. The well-known 
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Figure 2. Evolution of the reduced front profile. Fixed coupling: left; Running coupling: right. 
The reduced front profile (k 2 /Q 2 a ) lc M(k/Q S (Y),Y) is plotted against log(fc 2 /Qs) for different 
rapidities. The various lines correspond to rapidities from 2 (lower curves, full line) up to 10 
(upper curves). Note the similarity of the wave fronts, but the quicker time evolution (in s/t) for 
fixed coupling, by contrast with the slow time evolution (in t 1 / 3 ) for the running coupling case. 

geometric scaling property |10j is obtained for the solution of the non-linear equation 
(|3*|) at large enough energy. In our notation, the geometric scaling property can be 
written 

N{Y,x m ) =N(x 01 Q s (Y)) , (6) 

where 

Ql(Y) = exp L^Y - A logy - 3 / 2*1 + 1 
{ 7c 2 7c ( 7c ) 2 y ax" (7c) VY J 

plays the role of the saturation scale squared. Note that the solution (J5J mathemat- 
ically requires an value initial condition A/"o(fc— >+oo, Yq) <C l/k 2lc which is realized 
in first order QCD by color transparency TVo ~ 1/k 2 . Note that the result gives 
a rigorous proof of previous evaluations based on linear evolution with boundary 
conditions (first term: [TT], second term: the third is new [7]. 

It is possible to show that the result (J5J is more general, by various extensions 
of the F-KPP solutions. First, some general arguments confirmed by numerical 
simulations (see the review in 9 ) lead to expect the same result for the full nonlinear 
equation PJl. It is independent of the precise form of the non- linear damping terms 
and from the initial conditions (provided the transparency condition is fulfilled). 
Hence the "Universality" property. Second, the results can be extended to running 
a. One interesting difference with the fixed a case is the "abnormal" diffusion 
approach to scaling (in t 1 / 3, ~ y 1 / 6 instead of i 1 / 2 ~ Y 1 / 2 ), see Fig.©. 
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